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Abstract 

A new "bond-algebraic" approach to duality transformations provides a very 
powerful technique to analyze elementary excitations in the classical two- 
dimensional XY and p-clock models. By combining duality and Peierls ar- 
guments, we establish the existence of non-Abelian symmetries, the phase 
structure, and transitions of these models, unveil the nature of their topo- 
logical excitations, and explicitly show that a continuous U(l) symmetry 
emerges when p > 5. This latter symmetry is associated with the appear- 
ance of discrete vortices and Berezinskii-Kosterlitz-Thouless-type transitions. 
We derive a correlation inequality to prove that the intermediate phase, ap- 
pearing for p > 5, is critical (massless) with decaying power-law correlations. 
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1. Introduction 

In this article we investigate, via the use of dualities, two-dimensional 
(D = 2) classical systems, such as the XY and clock models [1, 2], that 
display Berezinskii-Kosterlitz-Thouless (BKT)-type transitions [3, 4, 5, 6]. 
BKT transitions, notably characterized by essential singularities in the free 
energy, emerge in many physical situations including screening in Coulomb 
gases [7], surface roughening [8], melting in D = 2 solids [4, 9], and many 
other classical and quantum problems, such as deconfinement in D = 3 + 1 
lattice gauge theories [10, 11, 12, 13]. We study these models by invoking 
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a bond- algebraic approach that we have recently developed [12, 13, 14, 15]. 
Within our duality-based method, one relates singularities in the free energy 
at one temperature (or coupling constants) to those at a dual temperature 
(or dual coupling constants). 

Specifically, we investigate exact dualities of the D = 2 XY and p-clock 
models, and exploit those dualities to unravel their phase structures. Those 
transformations are exact even for finite systems after appropriate boundary 
terms are included. It is noteworthy that unlike nearly all of the analytical 
work done to date, our dualities do not rely on the approximation scheme 
of Villain [16, 10, 17, 18], yet they can be related to the exact dualities of 
the Villain model in appropriate limits. Furthermore, our analysis leads to 
exact dualities for general p-clock models and yields a better understanding 
of the appearance of two transitions in systems with p > 5 states (the XY 
model with only one transition is recovered in the p — > oo limit). By fusing 
our duality results with the Peierls argument, we will be able to (1) prove 
that p > 5 clock systems can be made to be self-dual; (2) prove by a Peierls 
argument that there exists a lower ordering temperature ~ 1/p 2 , asso- 
ciated to domain-wall excitations, below which the global Z p symmetry is 
broken; (3) demonstrate that a second transition occurs at a temperature 
~ 0(1) when p > 5 (in the self-dual case, it follows that if is not 
the self-dual temperature T* , then there must necessarily be a second phase 
transition at of an identical character); and (4) characterize the nature 
of the topological excitations, and further explain that at p = 5 a new type of 
topological excitation, with an associated discrete winding number, appears. 
Our considerations suggest that these discrete vortices may proliferate above 
the temperature T^ 2 \ We will also (5) determine an analytic expression for 
the self-dual temperature T*, which relates and clarifies temperature scales 
discussed in [19], and most importantly, (6) establish the non-Abelian poly- 
hedral symmetry group P(2,2,p) of the p-clock and related models, and 
explicitly unveil the U(l) continuous symmetry that emerges when p > 5. 
Indeed, the latter is intimately tied to the existence of the BKT transition. 
Finally, (7) we derive a correlation inequality to prove that the intermediate 
phase, appearing for p > 5, is critical (massless) with decaying power-law 
correlations. 

Despite several analytic [10, 20, 21] and numerical calculations [22, 19, 23], 
the precise nature of the two phase transitions (p > 5) is not completely un- 
derstood. It was proven [21] that for large enough p, clock models exhibit 
a BKT-type transition (actually, it has only been proved that there exists 
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an intermediate critical phase with power-law correlations). The question 
whether that still holds when p = 5 remains open [22]. By relying on ex- 
act results, we shed light on the character of the BKT transitions in these 
systems. We will relate the BKT transition to a continuous U(l) emergent 
symmetry of an usual type. Although BKT transitions are often discussed 
in terms of specific anomalous exponents and jumps in the helicity modulus, 
we will not address such non-universal issues. 

Our treatment of classical dualities is based on a new approach developed 
in Refs. [14, 13] that relies on the transfer matrix or operator formalism 
[11]. In statistical mechanics, two models a and b are dual if their partition 
functions Z a = tr [T^] and Z b = tr [T b N ] are related as (N is the linear size 
of the system in D = 2) 

Z a [K\=A(K,K*)Z h [K*], (1) 

with A some analytic function of the set of couplings K of model a, and 
dual couplings K*. In principle, Eq. (1) establishes an extremely broad 
relationship that could be achieved through many transformation schemes, 
including the standard one based on taking the Fourier transform of in- 
dividual Boltzmann weights [17, 24]. However, it was discovered in Refs. 
[14, 13] that low-temperature (strong coupling) /high-temperature (weak cou- 
pling) dualities correspond to a unitary equivalence of transfer matrices or 
operators, T a and T b , 

T b = U d T a Ul (2) 

with Wd a unitary operator. This observation is extremely insightful because 
there is a simple and systematic way to look for unitary equivalences be- 
tween physical operators, based on the notion of bond algebra, or algebra of 
interactions [14, 13]. 

The outline of this article is as follows. In Section 2 we define the classical 
XY model, and then in 2.1 establish its transfer operator. In Section 2.2 
we discuss the form of the exact one-dimensional quantum analogue whose 
partition function is that of the D = 2 classical XY model with coupling 
constants K\ and K 2 . In the limit of large coupling K 2 along columns, 
the quantum model is the 0(2) quantum rotor model. In Section 2.3, we 
establish the duality of the D = 2 XY model to a solid-on-solid-like and also 
to a lattice Coulomb gas-like models and, moreover, determine the disorder 
variables. These dualities do not rely on the Villain approximation scheme 
but are exact dualities obtained by our bond-algebraic method [13]. 



3 



We next proceed to analyze in Section 3 the p-clock model [1]. This model 
provides a particular controlled limit to the XY model (the p -> oo limit). 
We replicate the same steps undertaken in the analysis of the XY model 
of Section 2, but now the Weyl algebra [25] and the theory of circulant 
matrices [26] play a key role. We construct in Section 3.1 its transfer matrix, 
and proceed in 3.2 to establish the corresponding one-dimensional quantum 
Hamiltonian, that is not self-dual for p > 5. We study the dualities of these 
systems in Section 3.3. The system is exactly self-dual for p = 2,3,4, and 
becomes approximately self-dual for K2 K\ when p > 5. In Section 3.4 
we introduce a variant of the classical p-clock model that is exactly self-dual 
for all p. We examine, in Section 3.5, the exact and emergent symmetries 
of these systems and, notably, unveil the U(l) symmetry that emerges when 
p > 5. This continuous emergent symmetry is responsible for the existence 
of the intermediate critical (massless) phase. 

Finally, in Section 4 we utilize our previous findings to better understand 
the phase diagram of the p-clock model. Here we present an analytic expres- 
sion for the self-dual temperature T*, an important scale in the problem, and 
advance a Peierls argument. We also introduce a topological invariant, that 
we call the discrete winding number k, to unravel the nature of the topologi- 
cal excitations. Starting at p > 5 a new type of topological excitation appears 
with a non-zero value of k that one may call discrete vortex, and which is 
responsible for the phase transition to a disordered state. Domain-wall topo- 
logical excitations are key at low-temperatures and their energy cost depends 
on p, and on the relative spin configurations (except for 2 < p < 4). By using 
both duality and energy-versus-entropy balance considerations we show that 
the transition from the critical to the disordered phase scales as ~ 0(1), 
while the one from the broken Z p symmetry to the critical phase goes as 
~ 1/p 2 for large p. We derive a correlation inequality, i.e., show that 
the two-point correlation function G is a monotonically decreasing function 
of temperature, allowing us to prove that, for p > 5, the intermediate phase 
is critical (massless) with decaying power-law correlations. The appendices 
provide technical developments including a duality of the XY model to q- 
deformed bosons, illustrating the key physical difference between compact 
and non-compact degrees of freedom. 
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2. The XY model: A paradigm of BKT phenomenology 

The D = 2 classical XY model is the paradigmatic example of a system 
displaying a BKT transition at a finite temperature Tg^ T > 0. This model, 
also known as planar rotator or planar 0(2), consists of an N x N array of 
classical two-component spins S r located at the vertices r = ie\ + j e 2 
being integers) of a square lattice with unit vectors e M , \x = 1, 2, as indicated 
in Fig. 1. Its partition function is 

Zxy[Kn,h] = ^ ex P 

{Sr} 

where the spin S r = e\ + S% e 2 , the coupling = fij^ is the product of 
the inverse temperature (3 = l/k B T and the exchange coupling J M , and h is 
a temperature-rescaled external magnetic field. Fixing the magnitude of the 
spin variable S% to one allows us to re-write the partition function as 

ZxxlK^h] = ^ ex P 

where the continuous angle variables 9 r take values in the interval 9 r = 9ij G 
[0, 27r), i.e., it is a compact variable. We assume, without loss of generality, 
that h — hei. The sum over configurations represents an integral 




In the remainder of this article, we will concentrate on the case with zero 
external magnetic field, i.e., h = 0. The XY model displays a (global) con- 
tinuous U(l) symmetry, which amounts to the invariance of the model under 
a simultaneous rotation of every spin in the lattice by the same angle. In 
low dimensional systems with continuous symmetries, such as the XY model 
above, long-range order is more fragile. Thermal fluctuations may induce in- 
stabilities with the end result that long-range order is actually non-existent 
in D < 2 dimensions. Spin-wave excitations are responsible for destroying 
such an order. The Mermin- Wagner theorem formalizes this qualitative pic- 
ture. In the context of the D = 2 XY model, the theorem states that this 
system does not display spontaneous magnetization at finite temperatures. 
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Figure 1: The two-dimensional classical XY model. On each vertex r — i e± + j e-2, of the 
square lattice there is a classical spin S r = S% e\ + Slj. e2 of magnitude 1, i.e., S r ■ S r = 1, 
and = cos^, 5^ = sin^ with 9 r G [0,27r). Nearest neighbor spins interact with an 
exchange constant value J\ or J2 depending on the spatial direction. 



A common physical mechanism behind the formal proofs for both classical 
and quantum versions of the XY model, can be found in Ref. [11]. 

A phase transition is said to occur whenever a thermodynamic function 
of the system under study displays a non-analyticity. The latter may occur 
even when the ground state is unique, so that there is no spontaneous sym- 
metry breakdown. This is the case for the D = 2 XY model, that is known to 
have a special phase transition at a finite, non-zero temperature T^ T . This 
BKT transition is characterized by an essential singularity in the free energy 
and correlation length at Tg^ T . If T > T^ T , the correlators of the XY model 

decay exponentially with distance, as is typical of a disordered, paramagnetic 

(2) 

phase. In the low-temperature phase, T < Tg KT , the correlators decay alge- 

(2) 

braically with distance, just as if every temperature below T BK ; T represented 
an ordinary critical point. The fact that this power- law behaviour extends 
over the finite temperature range < T < Tg^ T , without long-range order, 
is known as quasi-long-range order. 
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2.1. A transfer operator for the XY model 

In this section we set up a transfer operator for the XY model, in prepa- 
ration for the detailed study of its duality properties and symmetries. We 
assume open boundary conditions in the ex-direction and periodic ones in 
the e 2 -direction. 

Since we are considering the XY model on a square lattice of size N x N 
we will need the operators 

L z>i =-\-^-, and e ±i§ \ i = 1, 2, • • • , iV, (6) 
satisfying the following commutation relations 

[L, ti , e ±i *']=±6 id e k,$ i. (7) 
The eigenstates of the unitary operators e ±id \ 

fF 6i \6 i ) = e* : * i \0 i ), ^g[0,2tt), (8) 

satisfy 

(e' i \6 i ) = 5{6' i -e i ), [* do, \e t )(e t \ = t. (9) 

J 

The plane wave eigenstates \rii) of L z ^ form an orthonormal basis of the 
Hilbert space of square integrable functions on the cirle £ 2 (U(1)) = Hi, and 
are related to the states \6i) via (9i\rii) = e' 9int / y/2ir. 

On the one hand, e ±l6 ' 4 represent position operators for the spin at site i, 
since their simultaneous eigenstates \6j) specify one, and only one point on 
the unit circle. On the other hand, L z ti represents their canonically conjugate 
momentum, the infinitesimal generator of translations 

e - iSL *.%) = \9i + S). (10) 

This last equation follows from Eq. (7), since 

The product states, 

\o) = <g>m, (12) 
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that are simultaneous eigenstates of all the position operators e ±l6, % form an 
orthonormal basis of the total Hilbert space Ji = 7^^. 

We have now all the ingredients needed to write down a transfer operator 
for the XY model. Consider for concreteness the following row-to-row (j to 
j + 1) matrix elements of the transfer operator 



(0'|T 2 |0) = exp 
and the diagonal operator 



A 1 



2 COSl^jj+i 

i=l 



(13) 



T 1 \6)=exp 



'N-l 



^2 K i cos(0 i+ ij - 9 it 



i=l 



\0). 



(14) 



Both matrices are defined in the basis of states introduced in Eq. (12). It is 
straightforward to check that if we set Txy = T 2 Ti, then 



tr =Z X y[K ll ,h = 0], 



(15) 



recovers the partition function for the XY model of Eq. (3), provided that 
we set the external magnetic field h to zero. 

Next, we rewrite the operators T 1; T 2 in terms of the operators introduced 
in Eq. (7). The result reads 



N-l 



N 

n 



2tt 



d9 e 



K2 cos I 



(16) 



as can be checked by taking matrix elements of T 2 Ti in the basis of Eq. (12). 
Notice that T\ factors into a product of two-body operators that involves 
only nearest neighbors, while T 2 factors into a product of one-body operators. 
This important simplification is a direct reflection of locality. 

The relevant symmetries of the classical XY model translate into unitary 
transformations that commute with Txy = T 2 Ti. Besides the obvious ge- 
ometrical symmetries of the lattice, Txy commute with two operators that 
represent internal, global symmetries. The continuous global U(l) symme- 
try under global rotations of the classical spin direction 9 r — > 9 r + a, W, 



guarantees that [L z , Txy] = 0, where L z = ^i=i -^-M ^ s ^ e total angular 
momentum. There is also a discrete symmetry Cq = Yli=i Coi that is alluded 



S 



to, somewhat inaccurately, as "charge conjugation" [27]. The operator 
acts on position eigenstates as C^Oj) = |27r — 0j), and on angular momentum 
eigenstates as C^rii) — \ — rii). Thus, = 1 and = Coi. Since Co does 
not commute but rather anticommutes with L z , we notice that the full group 
of internal symmetries of the XY model is non-Abelian. 

2.2. Hamiltonian form of the XY model 

Models that can be written in terms of an Hermitian transfer matrix or 
operator, such as the XY model, can be translated into quantum-mechanical 
problems [11] by simply defining a quantum Hamiltonian according to 

Hxy = — ln(Txy), or equivalently, Txy = e~ Hxy . (17) 

While this is a powerful tool, often used in the literature, its actual value 
is diminished by the technical problem of computing ln(Txy) and, perhaps 
more importantly, because Hxy turns out to be a highly non-local opera- 
tor. The standard way out of this difficulty is to make approximations that 
solve both of these problems. The qualitative picture that emerges is intu- 
itively appealing but problematic if the approximations are not reasonably 
controlled. 

We will not determine Hxy in closed form, but rather we will compute 

H^ = -\nT^ M = l,2. (18) 

in closed form, and then exploit the Baker- Campbell-Haussdorf (BCH) for- 
mula to obtain the expansion 

H = - \n(e- H2 e- Hl ) = H 1 + H 2 + [H 1} H 2 }/2 + ■■■ . (19) 

Since Txy — T 2 T 1; as defined in the previous section, is not Hermitian, we 
have to set Hxy = {H + H j; )/2. This only affects terms quadratic and higher 
order in the commutators. We will also study the conditions for the non- 
diagonal (kinetic) part of the Hamiltonian H 2 to reduce to the intuitively 
appealing form H 2 oc ^ L 2 zi /2. 

Referring back to the operator form of T 1; T 2 , Eq. (16), we see that it 
is straightforward to compute H\. Since Tx is already in diagonal form, we 
obtain 

N-l 

#i = -X^i C0S (4+i-^). (20) 

i=i 
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On the other hand, T 2 is not diagonal. Computing H 2 is further simplified by 
the fact that T 2 factors into the product of N commuting one-body operators. 
It follows that H 2 = ; w ith 

e -H 2li / dQ e K 2 cos9 e -WL z ^ ^ 

Jo 

Next we notice that, since e~' 9lLz ' 1 e~' 6,2Xz ' i = e ~'( 6 ' 1+£ ' 2 ) Lz ' 1 , H 2ji should be of 
the form (see Appendix A) 

H 2ji = - dQ a K2 {6)e- WL >>\ (22) 
J 

By combining this expression with Eq. (21), we get an equality between 
functions of L Z)i that can be evaluated in that operator's diagonal basis. The 
result is an infinite set of equations 

/ dO a K2 {6)e- Wn = \n(2irI n (K 2 )), n e Z, (23) 
Jo 

that one can use to determine aK 2 {9) by taking the Fourier transform of the 
left-hand side. The modified Bessel functions of the first kind and of integer 
order n, I n (K 2 ), satisfy [28] 

00 

e K 2COS e = Y^I n (K 2 )e< en = I (K 2 )+2^2l n (K 2 )co8(6n), (24) 

neZ 71=1 

after using the relation I_ n (K) = I n (K) [28]. 

It follows from Eq. (23) that ct#; 2 (27r — 9) = ax 2 (#) ; so that we can make 
the substitution e *>* — > cos(#L 2) j) in Eq. (22). Then, we can Taylor- 
expand the cosine function to get 

00 

H 2 , = -J2 a m(K 2 )L%, (25) 



m=0 



with a m (K 2 ) = (-l) m J^d9 6 2m a K2 (6) / '(2m)!. This equation provides a 
very convenient representation of H 2 ^, especially if we are allowed to discard 
terms beyond m — 1. To investigate this possibility, consider Eqs. (21) and 
(25), and evaluate those expressions in the L z /s diagonal basis to get 

00 

a m (K 2 ) n 2m = ln(2vr/ n (^ 2 )). (26) 

m=0 
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For large K 2 , the functions I n (K 2 ) have the following asymptotic expansion 
[28] 

^) ~ £ ( " ir ^# (27) 



m=0 



2 



with 



, N , , , (4n 2 -l)(4n 2 -3 2 )---(4n 2 -(2m-l) 2 ) 

c n = 1, c m n = ^ ^ ^— ^ ^ U. m>\. 

ml 8 m 

(28) 

Notice that this can be rearranged into an expansion in n 2 that can be 
compared to the left-hand side of Eq. (26). Keeping, for each m, only the 
leading order in I/K2, and expanding the logarithm accordingly (ln(l + x) ~ 
x), we obtain 

ln(2nI n (K 2 )) ~ In ( J — e K A + £ (-l) m 0m „,^ m . (29) 

' m=l 



V V 7 ^ 2 m m!Js: 9 

Comparing with Eq. (26), 

/ / 2n \ (— l) m 

flo(iG)~ln^— e^J, M^)~^^, m>l, (30) 

so that a m+ i/a m ~ — l/(2(m + 1)^2) ■ in summary, in the large limit, 

JV . AT-l 

iJxv ~ -iVa (^ 2 ) + ]T 2K~^* ,i ~ S Kl cos ^ +1 ~ ( 31 ) 

i=l 2 i=l 

which is the one- dimensional 0(2) quantum rotor model. 

5. Duality of the XY model without the Villain approximation 

We now exploit the detailed understanding we have gained on the ex- 
act operator structure of the XY model to look for its dual representations. 
The standard approach to the dualities of the XY model starts by replac- 
ing it with the Villain model (see Appendix C), then mapping the latter 
to the solid-on-solid (SoS) model, and finally mapping the SoS model to a 
lattice Coulomb gas [11]. In contrast, in this section we establish directly 
exact dual representations of the XY model. In Appendix B, we establish 
a duality to a g-deformed boson Hamiltonian which illustrates the fact that 
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non-canonical bosons need to emerge because of the compact nature of the 
degrees of freedom of the XY model. 

Our methodology starts with the transfer operators 7\, T2 introduced in 
Eq. (16). The algebra of interactions, or bond algebra in the language of 
Refs. [12, 13, 14, 15], underlies their basic structure. In this case, this is the 
von Neumann algebra Axy generated by the bonds 

L zi , e^-^\ i = 2,---,N. 



Notice that T±, T 2 G -4xy, since these operators are expressible as sums 
of products of the bonds listed in Eq. (32). Axy reflects the interactions 
present in the XY model and is at the same time easy to characterize in 
terms of relations. Then we can look for other dual realizations A^ Y that 
are isomorphic images of Axy- By the general properties of von Neumann 
algebras, it must be that A^ y = U^AxyU\, with U.& unitary, provided both 
algebras act on state Hilbert spaces of the same dimensionality. This is all 
we need to establish a duality for the XY model. The dual partition function 
is determined from the dual transfer operator = U^TxyU.\. 

The goal of this section is to look for a dual representation of the XY 
model in terms of integer-valued degrees of freedom, so we can expect the 
dual bond algebra A^y to act on the state space ®f =1 £ 2 (Z). Let us introduce 
the states |n) = <S) i=1 \rii), and the operators Xi, Ri 

Xi\n) = rii\n), (32) 
Ri\n) = I • • • ,n;_i,n; - l,n i+1 , • • •), R\\n) = \ ■ ■ ■ ,ni_i,rij + l,n i+1 , • • •), 

that satisfy the algebra 

[Xi, Rj] = -SijRj, [X h R]} = 5 ij3 R], RjR] = 1. (33) 
Then, the operators 

X u X i+1 -X u R u R\, i = l,.-. ,N-1, (34) 

generate an isomorphic dual representation ^.^y °f ^ ne bond algebra of the 
XY model. The isomorphism $d connecting the two bond algebras 

L z ,x — ^ Xi, L z i — h Xi — i — 2, • • • ,N, (35) 

e+ Kk + ,-k) 1^ Rh e-MJ^Rt z = l,-..,JV-l, 
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L z ,l e +\{6 2 -6i) L z-2 e +Kk-9 2 ) L Z;3 



X — x — X 




Figure 2: The bond algebra isomorphism $d(-) = U&-U A defined in Eq. (35), illustrated 
for N = 3 sites. 



is illustrated in Fig. 2, for N = 3 sites. The resulting dual transfer operators 
are 

N-l 

T± = JJ e ^+^), (36) 



i=l 
2tt 



T. 



D 



d9e 



K 2 cos6-\X 







N , 2 tt 

1 n / dd 



i=2 







The next and last step is to compute the dual partition function 



2% = tr [(TfT?)*] = tr [(75° > 



in the product basis states of Eq. (32). 

T 2 D is already diagonal in that basis, leading simply to 



T 2 D \n) = exp 



N 



where 



V K (n) = - In / cifle 



2- 



,if cos Sg— i9n 



ln(27r/„(AT)). 



(37) 



\n), (38) 



(39) 



The evaluation of (n'\Tf\n) factors into the computation of the one-body 
matrix elements, 

(vf^^+^lm). (40) 

This is simplified by noticing that the Fourier transform operator of Eq. 
(A. 2) maps 



F} RjFi = e 



F}B\F t = e" 



(41) 
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thus putting the exponential in diagonal form. It follows that 

(n^e^^+^Vi) = 27re- VK ^+ 1 - n ^ ) . (42) 

If we now put all the pieces together, we arrive at the conclusion that we 
have obtained the exact duality 

Z Xvl K n] ST p [-Ef =1 E^ 1 (^ 2 K+i, J -n l ,,)+y Xl (n l , J+1 -n l , J ))+Ef =1 VK 2 (ni, J )] 

(43) 

This illustrates a typical characteristic of dualities: the coupling K\ (K 2 ) 
in the ei(e2)-direction of the XY model regulates the interaction in the 
orthogonal e2(ei)-direction of the dual model. 

It is standard to argue that the Villain model is an excellent approxima- 
tion to the XY model, specially at low temperatures. As shown in Appendix 
C, the Villain model is dual to the SoS model. Thus, it must be that the SoS 
model is approximatly related to Z^ Y defined by the right-hand side of Eq. 
(43), at least for low temperatures. Consider then the limit of large K%, K 2 
(i.e., low temperatures). We can then use the asymptotic expansion of Eq. 
(27) to obtain an asymptotic form of the dual potential of Eq. (39), 

V K (n)^ — + c(K), K^oo, (44) 

where c(K) is independent of n and can be computed from Eqs. (39) and 
(27). It follows that 

XY / j 

{ n i,j} 

(45) 

to the same level of approximation. Thus we have recovered the well known 
result that the XY model at very low temperatures (strong coupling) is well 
represented by the (approximately dual) SoS model at very high tempera- 
tures (weak coupling). 

The action of the duality of Eq. (35) can be extended to act on the 
operator e' &N as $d(e~' e ' JV ) = Rn- It follows that <3>d generates the following 
set of N dual variables (below we distinguish a dual variable by an overtilde), 

. N 

e -A = $ d ( e -i**) = JJ R m) ; = !,... ,jV (46) 



L z ,i = $d(L z>1 ) = X h L z>i = $ d {L z ,i) =Xi- i = 2, • • • , N. 
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The dual variables satisfy the algebra of Eq. (7), confirming that $d de- 
fines an algebra isomorphism. Since this is also the algebra of the variables 

Ri, rI, Xi, we see that the dual variables e ±] ®\ L 2 i afford an alternative rep- 
resentation of the elementary degrees of freedom. But what is their thermal 
behavior? This crucial question can be answered easily because $d amounts 
to a unitary transformation. It follows that 







11 I j xy e i xy e i XY ) 

tr (T^ N ~ n ~ s) e ,6 ™+r T X y e|&» T£ Y ) 
~ZP 



that should be compared to 



um-k-r.tl+s 



1 V XY ^m+r J-XY -^m ^XYJ 



■Wm,») j (47) 



/ \ _ - V^XY — m+r -XY — m -XY/ /.on 



z 1 

^XY 

The classical dual variables e _l6 ' m '" are difficult to define directly, but they 
are well defined in the sense that any correlator 

in the ensemble Z XY can be computed by a straightforward generalization of 
(47). 

The duality to a lattice Coulomb gas is of a very special nature (Poisson 
duality) and quite different from every other duality discussed in this article 
(or in literature on dualities in general). Its general features are discussed 
in Ref. [13]. Here, we briefly summarize the exact Coulomb gas-like dual 
model for the exact dual partition function Z XY computed above. According 
to Ref. [13], the lattice gas representation of Z XY is defined by 

^XY D = £ e ^ {M ' (50) 

{n r } 

with interaction energy £ D defined through the global Fourier transform 

e-K 1 = Jn tre =,E.» e E.E,«v„ ( , W( ,,, (51) 
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The interaction energy 8 can be determined in closed-form in the limit in 
which Z£y reduces to the SoS model, and corresponds to the standard lattice 
Coulomb gas result [11]. We see from Eq. (51) that Poisson dualities are in 
general only of practical use for models with Gaussian energy functionals. 



3. The p-clock model: A close relative of XY 

The p-clock model, also known as the vector Potts or Z p model, represents 
a hierarchy of approximations to the XY model as a function of the positive 
integer p, and is a test ground for rich critical behavior. In D = 2 dimensions, 
the configurations of the classical p-clock model are described by a set of p 
discretized angles 9 r 

6 r = s r = 0,l,--- ,p-l. 52 

p 

The partition function is given by 



z p\ k A = Y exp 

{0 T } 



(53) 



This is, in appearance, identical to the classical XY model (h — 0), except 
for the essential fact that the degrees of freedom {6 r } are now discrete and 
countable. In the large p limit (p — > oo), however, one supposedly recovers 
the XY model. Since the p points e' dr close a Z p subgroup of U(l) (the group 
of p'th roots of unity), the p-clock model manages to provide an approxima- 
tion to the XY model that features a finite number of states per site, without 
sacrificing the XY's natural group structure. Also like the XY model, the 
p-clock has interesting but hard to uncover duality properties [20]. We will 
address this problem by the same methods applied to the XY model. In 
fact, we will follow the discussion of the XY model as closely as possible, to 
highlight the connections between the two models. 

3.1. A transfer matrix for the p-clock model 

To introduce a transfer matrix for the p-clock model, we need to define 
a suitable Hilbert space and a set of basic kinematical operators. Let Z p be 
defined on an N x N square lattice with open boundary conditions on the 
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ei-direction and periodic ones on the e 2 -direction. The states on each site 
% — 1, • • • , N of a row can be described by orthonormal vectors 

\st), Si = 0, ••• ,p-l, (54) 

such that Si represents the discrete angle 2iiSi/p. They span the state space 
H p ,i at site i, so the total state space is just H p = (£)f =1 Ti p ,i- If we write 
\s) = (£) i=1 | Si) for elements of the product basis of H p , then we can define a 
matrix 

(s'\T 2 \s) = exp 

related to any pair of adjacent rows j, j + 1, and the diagonal matrix 
Ti\s) = exp 

These definitions guarantee that Z P [K^] = tr [(Tp)^] = tr [(T 2 Ti) N }. 

The degrees of freedom of the p-clock model (at any one site of the lattice) 
can take any value out of a discrete, equidistant subset of points of the 
unit circle. To proceed in analogy to Section 2.1, we need to introduce 
position operators and their conjugate momenta in this discrete setting. The 
formalism that emerges was used extensively by Schwinger in his work on the 
foundations of quantum mechanics [25]. In what follows, we consider only 
one site (one degree of freedom), for the sake of clarity. We will consider all 
iV sites again near the end of the section. 

It is easy to restrict the position operators e ±id used for the XY model to 
the subset of configurations available to a clock handle in the p-clock model. 
The result is the operator U satisfying 

U\s) = oo s \s) , s = 0,---,p-l, (57) 

with to = e l27r / p representing a pth root of unity. The position operator U and 
its Hermitian-conjugate W satisfy UU' = 1 = U p . The momentum operator 
V conjugate to U rotates any state counter-clockwise to its nearest-neighbor 

V|0) = b-1), V|1) = |0), V \p-l) = \p-2). (58) 



N 



5> 



2 COS 



27r Si j+i 2nsij , 



i=l 



p 



(55) 



N-l 



E ,2-KSi+ij 2-KSu, 
Ki cos( 



i=l 



p 



p 



\s). 



(56) 



17 



Momentum and position operators are represented, in (p x p) matrix form, 

as 



V 



/o 1 o ••• o\ 

1 ••• 

•■■ 1 

••• 0/ 



and U = diag(l, u, 



(59) 



It follows that implements a clockwise rotation, and that VV' = 1 = V p . 
The fundamental algebraic relation 



VU = ujUV 



(60) 



follows directly from the definitions of U and V. 

As is well known from quantum mechanics, the ordinary position operator 
x and its conjugate momentum operator p are related by a Fourier transform 
J 7 , a unitary transformation in the space of wave functions. Essentially the 
same holds for the operators U, IP and V, V'. The appropriate unitary 
transformation in this context is the discrete Fourier transform F, that in 
matrix form reads 



F* = ± 



(1 


1 


1 


1 




u 2 


1 


u 2 





CO 



1 

p-1 



W 2(P-1) 



(61) 



\1 LO^ 1 LO^ 2 ■■■ 

This is also known as Schur matrix [29]. It follows that [25] 

FUF ] = V\ FVF ] = U, 

and so the eigenvectors of V, s = 0, 1, • • • ,p — 1, 

p-i 



a; s 



1 

) , with \s) = — ^u ss \s), 



(62) 



(63) 



are easily determined via a Fourier transform of the eigenvectors of U. 
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In the mathematical literature V is known as the fundamental circulant 
matrix. This is so as it generates the algebra of circulant matrices [26] (mean- 
ing that any circulant matrix C is of the form C = Y^n=o a mV m , a m G C). 
Together, U and V generate the full algebra of (p x p) complex matrices [25], 
that we continue to call the Weyl group algebra, to emphasize that we are 
working with a distinguished set of generators. This shows that they consti- 
tute a convenient basis set of kinematic operators, because we can write any 
other operator in terms of them. 

We need to re-introduce the row spatial index i to apply the technology 
just developed and rewrite the transfer matrices of Eqs. (56) and (55) in 
operator form. In what follows, Ui, LfJ, V$, V i , for % — 1, • • • , N, will be our 
basic set of operators. They commute at different sites, satisfy the relation 
(60) at any one site i, and act on the state space H p = (^)j =1 H p ,i- One then 
obtains 

N-1 N p-l 

Tl = II e^r {u >+i Ul+u * +lU >\ T 2 = Y[Y1 e I<2COs{2wm/p) V? m . (64) 

i=l i=l m=0 

This last expression for T 2 follows from the fact that (s^|V^ |s,) = unless 
s'i — Si = m modulo p (mod(p)). It should be compared to the analogous 
expression for the continuum circle, Eq. (16). 

3.2. Hamiltonian form of the p-clock model 

In this section we compute the Hamiltonian form of the p-clock model fol- 
lowing the strategy of Section 2.2. We start by computing = — lnT M , fi = 
1,2, with Ti, T 2 as defined in Eq. (64). 

Since T\ is diagonal, we can write 

N ~ l K 

Hi = ~Y, -^( U U U i + U}U»i). (65) 



2 

i=l 



H2 = YliLi ^2,1, on the other hand, is not as easy to write down. H2,i is 
defined as 

p-i 

e ~ H 2,i — g^2 cos(27rm/p) y^ m _ (Q6) 

m=0 

As explained in Appendix A, we can solve this equation to obtain 

p-i 

H 2ti = * m {K 2 )V? m , (67) 

m=0 
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with 

p— l „ / p— i 

a. 



^ s=0 ^ \Z=0 ^ / 

Then, the Hamiltonian if p for the p-clock model follows 

N-l K N p-1 

H P = ~Y1 -rWliUi + U}U i+1 ) EE M^O^, (69) 

i=l i=l m=0 

provided we truncate the BCH expansion of lnT p to linear order (see the 
discussion in Section 2.2). We notice for future reference that the discrete 
Fourier transform F = Yli=i Fi maps H p — > F^HF = H p , with 

N-l jyr N P-1 

H P = -Y, "y (KU^ + V?V i+1 ) EE a m{K 2 )Ut m . (70) 

i=l i=l m=0 

As discussed in Appendix A, the coefficients a m (K 2 ) have simple asymp- 
totic forms in the limit K 2 —> oo. The corresponding approximation to H p 
reads 

H p ^ Na (K 2 )-K 1 H u -2a 1 (K 2 )H v , (71) 

with 

N-l N 

Hu = -(U N + 4 + E^i+i^ + ' ^ = 2 E W + ( ?2 ) 

i=l i=l 

and (see Eq. (A.12)) 

ai(if 3 ) « e ^(coB(!hr/p)-i) j ao (j^ 2 ) ~ (73) 

Equation (71) shows a boundary term (—Ki(Un + U^)/2) not present in Eq. 
(69), and that we include to make this approximation to H p exactly self-dual 
[13]. 

The approximation made in going from Eq. (69) to Eq. (71), that keeps 
only V} and V^ p = Vi, is reminiscent of the one introduced in Section 
2.2 based on Eq. (26), whereby we replaced the operator cos(9L Zti ) for the 
simpler L 2 zi in Eq. (31). Indeed, the two aproximations coincide in the 
p — > oo limit. A simple way to see this is to notice that we can realize the 
operator V} directly in the Hilbert space of the XY model as V} — > e~' 2nL ~-^ p . 
Then, in the limit p — > oo, 

Vi + V? -> e a " L '- i/p + e -' ,2wL ** /p « 2 - (2n/p) 2 L 2 zl . (74) 
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m y] U 2 U[ y] U 3 U< 

Figure 3: The duality isomorphic mapping $d of Eq. (76), for N = 3 sites. 

3.3. Dualities of the p-clock model 

The dualities of the p-clock model appear as isomorphic representations 
of the bond algebras associated to the transfer matrices defined in Eq. (64). 
The bond algebra A p generated by 

Vi, Vl Vi, Vl UVU, V$Ui-i, i = 2,--.,N, (75) 

is simple to work with and adequate to our purposes. It has a dual (iso- 
morphic) representation A® generated by the same bonds listed in Eq. (75), 
except for Vn, V^ that have to be removed from the set of generators, and 
replaced by U±, U{. The duality isomorphism $d : A p — > A® reads 

Ul +1 U, t A Vl i = l,---,N-l, (76) 

V 1 -% U u Vi -% uUUi, i = 2,---,N, 

together with the corresponding Hermitian-conjugate entries (since it must 
happen that ® d (0^) = ® d (Oy). $ d is illustrated in Fig. 3, and should be 
compared with the duality of Section 2.3 for the XY model, illustrated in 
Fig. 2. 

The dual form Z® = tr [(T 2 D T 1 D ) Ar ] of the p-clock model that follows from 
Eq. (76) is defined by the dual transfer matrices 

JV-l 

T ? = U e ^ (Vi+VI) > ( 7? ) 

i=l 

p-1 N p—l 

rpD _ ^ e K 2 cos{2%m/p)jjfm x J~J ^ £ K 2 cos(2nm/p) (JjtjJ^^™ 

m=0 i=2 m=0 
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Clearly, the p-clock model is not self-dual for arbitrary p and arbitrary 
couplings. However, the model is approximately self-dual in the extreme 
anisotropic limit with K 2 3> K\, and it is exactly self-dual for p = 2,3,4, 
and any coupling. We study these aspects of the p-clock model in the next 
sections and in Appendix D. 

We mention, in closing, that there is a p-state model that approximates 
the p-clock in the same sense in which the Villain model approximates the 
XY model. This Z p Villain model [10] is exactly self-dual for any p, but is 
otherwise quite different from the self-dual p-clock model to be introduced 
next. 

3.4- Self-dual classical p- clock model 

In this section we introduce a classical model Z sdp that we call the self-dual 
p-clock. It is closely related to the p-clock model (identical for 2 < p < 4), 
yet it is exactly self-dual for any value of p, and has the distinct advantage 
over the Z p Villain model [10] that its transfer matrix is remarkably simple. 
To define the self-dual p-clock, we introduce the transfer matrices 

N-l N 
i=l i=l 

(78) 

Then, Z sdp [a , a u Ki] = tr [T^ p ], with T sdp = T 2 [K 2 \T 1 [ao,a 1 \, and where a 
and a,i are free parameters of the model to be determined, for instance, by 
the requirement that the approximation 

p-i 

e K 2 cos (2nm/p) y\m ^ e a + ai (V l + V l t ) 

m=0 

be as good as possible for arbitrary K 2 . 

Z s d p is self-dual due to the existence of a unitary transformation that 
maps 

T^Iu] T 2 [0,/^], T 2 [a , ai ) e^'TiK], (80) 

This fact results from a bond-algebraic analysis, but we omit the details 
which can be found in Ref. [13]. It follows from Eq. (80) that the self-dual 
line is specified by K\ = 2a\. The next issue then is to understand the 
structure of Z sdp in terms of classical variables. On one hand, it is clear that 
the interaction energy in the ei-direction is still of the form K\ cos(27r(s i+1 j — 
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s i,j) Ip)- O n ^ ne other hand, the interaction energy in the e 2 -direction, it(sjj- 
Sij+i), is determined by the relation 

p-i 

3 a +o 1 (y i +v; t ) 



^ e «M]/.t m = e ao+ai(V>+y/)_ (g-g 
m=0 

Then, from Eq. (A. 8), 

1 b/2] 9 

|p/2]^ P ' ^ ' 

([p/2] denotes the largest integer smaller than or equal to p/2). With the 
interaction potential u defined in this way, 

£sd P [a ,ai,^i] = e^J [u(a «- 8i ^ +l)+KlC08(!hr(ai+1 ^- s< ^ )/p)] (83) 

is exactly self-dual under the exchange 2a\ -B- i^i. 

One can see directly from Eq. (82) that u(p — m) = u(m). It follows that 

u[m) = 2_^K 2t r cos( ), (84) 

with couplings 

[P/2 



1 n 

. 1 v-^ / 27rrm . . 
K 2 ,r = a d rfi + t-t-t ) j cos( ) In 



1 V r 2?rmg 2ai cos(2Wp) 



m=0 

(85) 

determined by Eq. (82) and the orthogonality relation of Eq. (A. 13). The 
point to notice is that to make the p-clock self-dual, we need to add higher- 
order harmonics (terms cos(27rm(sjj — Sij+i)/p), with m = 2, • • • , [p/2]) to 
the basic cosine interaction. We show next that Z s d p becomes a very good 
approximation to the standard p-clock model in a suitable limit. 

A comparison of the self-dual p-clock to the standard p-clock model shows 
that the latter has an approximate self-duality for p > 5. As explained in 
Appendix A, in the limit in which K 2 is large, Eq. (79) becomes almost 
exact, with 

a ^K 2 , ai « e K ^T-V (86) 



23 



(see Eq. (A. 12)), so that 



Z P [K,} « e^tr (e** H ve K * H ") 



2 ai H v KiH^N 



e N2K Hx 



/ 2aiH Ue KiH v \X 

e N2{K2 - K '^Z p [K*} (87) 



(see Eq. (72)), with dual couplings 



Kl=2e K «"»f- 1 \ K* 2 = ln( f l/2) . 

cos — — 1 

p 

We emphasize that this approximate self-duality, in the extreme anisotropic 
limit, is valid for any value of p. We consider exact self-dualities for the par- 
ticular cases p = 2, 3, 4 in Appendix D. 

3.5. Exact and emergent symmetries of the p-clock model 

Non-Abelian, discrete symmetries. The representation of the transfer ma- 
trix T p = T 2 Tx, Eq. (64), is very convenient for understanding the internal, 
global symmetries of the p-clock model. It is apparent that the model has 
an Abelian Z p symmetry, but, as it turns out, its full group of symmetries is 
considerably larger and non-Abelian, provided p > 3. To prove this we will 
show that there are two Hermitian operators 

N N 

Co = Y[Coi, Ci = IJ <7 U , (89) 

i=l i=l 

that commute with T p and satisfy 

Cl = (% = {C<£if = t. (90) 

These relations show that, if p > 3, Cq and C\ generate a unitary represen- 
tation of the so called polyhedral group P(2,2,p) [30] of order 2p, and so 
the group of internal symmetries of the p-clock model is at least as big as 
this non-Abelian group. Notice that C$Ci = Q, known as Z p charge, gener- 
ates a Z p subgroup of P(2,2,p). This is the standard Abelian symmetry of 
the p-clock model that gets broken in the low-temperature ordered phase (see 
Section 4). It becomes a U(l) symmetry in the limit p —> oo, and corresponds 
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to the usual continuous symmetry of the classical XY model 6 r — > 6 r + a, 
with a an arbitrary real number. 

As Co and C\ are products of one-site operators, let us focus on a single 
site % for now. We define the operators Coi and Cu by specifying their action 
on the basis of Eq. (54), 

C i\si) = | - s^, Cu\si) = |1 - Si), Si = 0, ••• ,p- 1. (91) 

The arithmetic in these definitions is modular, mod(p). For example, if p = 5, 
then C 0i |0) = | - 0) = |0), C 0i \l) = | - 1) = |4), and so on. Keeping this 
in mind, one can check that C and C\ are Hermitian, C^ = C 0i , C\ { = Cu, 
and satisfy the relations listed in Eq. (90) for Co and C\. In particular, as 
CoiCu = Vi, {CoiCuY = 1. If p = 2, then Co = 1, and C\ generates the 
Abelian Z2 symmetry of the Ising model which is different from the non- 
Abelian P(2,2,2). 

A routine calculation shows the action of Cqi, Cu on the discrete position 
Ui and momentum Vi operators, 

CoViCo = V?, C 1 V i Ci = V}, (92) 

CoVfio = U}, C 1 U l C 1 =cuUj. (93) 

It is easy to check that Co and C\ commute with T p . The operator Co is known 
in the literature as the "charge-conjugation" operator [27]. However, as we 
alluded to earlier, this is something of a misnomer. Geometrically speaking, 
Co is the exact analogue of the parity operator V\x) = \ — x) on the real line. 
In fact, Co is related to the discrete Fourier transform as F 2 = (i 7 ^) 2 = Co, 
just as its counterpart on the real line J 7 is connected to the parity operator 
as J= 2 = V. 

We wish to emphasize that these non- Abelian symmetries are shared by 
a large number of classical and quantum p-state models besides the p-clock, 
including the self-dual p-clock introduced in Section 3.4 and the Z p Villain 
models. 

Emergent U(l) symmetry. For p > 5 the discrete charge symmetry Q 
gets enhanced into a continuous U(l) symmetry. In reality, this is not an 
exact symmetry it is an emergent one [31], but it is essential to establish the 
intermediate BKT critical (massless) phase (see Section 4). Let us derive 
this emergent symmetry. 

Given the generators of the SU(2) algebra in the spin S = (p — l)/2 
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and S = (S + Y , one may study the transformation properties of the Weyl's 
group generators U and V under the U(l) mapping 



Us 



(94) 



Since U = oo~ U 2n / P i it commutes with Us- The transformation of V requires 
some thinking: Let us rewrite V, Eq. (59), as the sum of two operators 



V = V + A , with A 



/0 





\1 



o\ 




V 



(95) 



i.e., the matrix that has only a 1 in the lower- left corner. Then, 



''*V , and Us AU\ = e i(p_1) *A. 



(96) 



We are interested in analyzing how the transfer matrices Ti and T 2 of Eq. 



(64) transform under Us = YliLi^,i- ^ is indeed easier to analyze the 
Fourier transform transfer matrices, T M = F^T^F 



N-l 



N p-1 



f 1 =Y[ e^ {v ^ Vl+Vl+lV ?\ f 2 = Y[ ^2 e X2COs(2 ™ /p) f// m . 



(97) 



i=l 



=1 m=0 
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Clearly, U<p commutes with T 2 but does not with fi unless = 2ir/p, which 
not surprisingly corresponds to the (Fourier transform) discrete Q symmetry. 
However, is an exact continuous symmetry of the modified transfer matrix 

N-l 
i=l 

and becomes the usual U(l) symmetry of the XY model when p — > oo. This 
emergent symmetry may allow for the construction of spin-wave excitations 
in the critical region. Note that in the original transfer matrix represen- 
tation T\,T2, the continuous emergent symmetry is represented by FU^F^. 
Moreover, it is an emergent symmetry of both p-clock and self-dual p-clock 
models. 



4. Phase diagram: From the p-clock to the XY model 

We are thus left with the task of establishing the phase diagram of the 
p-clock model, the nature of its phase transitions and excitations, and its 
behavior as p — > oo. One may argue that the phase structure of the model 
is well understood [10] (see Fig. 4). At very low temperatures, there is a 
ferromagnetic phase characterized by long-range order of the two-point, spin- 
spin, correlation function G(\r — r'\) = (cos(# r — 6 r /)), and the breakdown of 
the Z p symmetry Q. At very high temperatures, the system is in a disordered 
phase with G(\r — r'\) decaying as an exponential function of the distance. 
For 2 < p < 4, these two phases are separated by a continuous second-order 
phase transition of the Ising (p = 2,4) or Potts (p = 3) type. (It is very 
easy to prove that the p = 4 case is identical to two uncoupled p = 2 Ising 
models [32], see Appendix D.) For p > 5 there is an additional intermediate 
critical phase separating the ferromagnetic from the disordered phase. It is 
characterized by a power- law behavior of G(\r — r'\) ~ \r — r'\~ v with a non- 
universal exponent rj, and by the absence of symmetry breakdown and quasi- 
long-range order. In the p — > oo limit the broken-symmetry phase disappears 
as one recovers the D = 2 XY model with a continous U(l) symmetry. This 
qualitative picture leaves several issues unresolved that numerical simulations 
have not been able to resolve either: 

• What is the nature of the two phase transitions for p > 5? 

• What is the nature of the relevant topological excitations in each phase? 
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What is the physical origin of the critical (massless) phase? 

What is the minimum p after which the transitions are of the BKT 

type? 



Zr, broken Disordered 
1 2 < p < 4 



2 nd order 



Zp broken Critical Disordered 
— X lllllllllll X 5 < p < oo 

BKT BKT 



Critical Disordered 
llllllllllllllllllllllllllllllll X P — > oo (XY 

BKT > 

a^/Kr (or T) 



Figure 4: Phase diagram of the p-clock model. For p > 5 there are three phases, the broken 
Z p (low-temperature) phase disappearing in the limit p — > oo (XY limit). A transition is 
of BKT-typc whenever it is associated to an essential singularity of the free energy. The 
critical phase is characterized by power- law correlations, i.e., quasi-long-range order, with 
non-universal exponents. 

To understand qualitatively the nature of the phases of the model, con- 
sider the ground state of the self-dual quantum Hamiltonian H p defined in 
Eq. (71), in the large (low-temperature) and small (high-temperature) K\ 
limits {a^K-i) = 0). Let us start with the broken Z p symmetry, low temper- 
ature sector that corresponds to the line (Ki,ax(K2) = 0). Then, the p-fold 
degenerate subspace of ground states is trivial to describe in terms of the 
simultaneous eigenvectors of the Ui of Eq. (57), 

l^o) = I [ \ s i)> with same s for all i. (99) 
i=i 

The ground state energy is Eq = —K\N for periodic boundary conditions 
(i?o = —K\{N — 1) for open boundary conditions), and (^qI^o) = $rs [33]. 
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The fully disordered, high-temperature phase is defined by the sector (Ki 
0, ai(K 2 ))- The ground state (E = —2ai(K 2 )N) is unique and given by 



N 



i=i 



(in terms of the eigenstates of Vi, Eq. (63)), and satisfies C |$o) = +| < ^o)- It 
is difficult to obtain exact results for arbitrary couplings. There is, however, 
an interesting exact relation that holds at the self-dual line K\ = 2a\(K 2 ) = 
K* and follows from the fact that the self-duality unitary U d becomes a 
new symmetry of the problem on that line. It is clear from Eq. (71) that 
H p [K*\ = -K*(H V + H v ), and U A H V U\ = H v , U A E V U\ = H v [13]. Since 
[H p [K*],h(d] = 0, we can choose the energy eigenstates \^ n ), n — 0, 1, • ■ -, to 
be also eigenstates of U d , U d \^ n ) = e'^ n |\l/ n ). Then 

E n = (* n \H p [K*]\V n ) = -2K*(V n \H v \y n ) = -2K*(V n \H v \V n ). (101) 

For 2 < p < 4, the p-clock model is exactly self-dual and the transition 
from the ferromagnetic to the disordered phase happens at the self-dual point 
K\ = 2a\(K2). For p > 5, Eq. (69) shows that the p-clock model is no longer 
exactly self-dual, but the self-dual approximation of Eq. (71) or (78) allow 
us to establish the following self- dual equation for arbitrary p 

bi = e ^(cos^-i) = i dm£ P (ai) ; = P y e 2a lCOs{ f m) _ {m) 

From the self-dual condition K\ = 2a\(K 2 ) one can determine the self-dual 
temperature T*. The self-dual point is a point of non-analyticity of the free 
energy for 2 < p < 4, but for p > 5 it is analytic. Some values are indicated 
in Table 1, assuming isotropic couplings Ki = K 2 = J/(/cbT). It follows 
from very general considerations (see Section 8 of Ref. [13]) that the two 
critical points c\ and c 2 bounding the self-dual point when p > 5 are exactly 
related by 



2a x {K 2 )\ Cl 2a x {K 2 



1. (103) 



It is interesting to analyze the large-p limit of the self-dual Eq. (102). In 
that limit 



lim -B p ( ai ) = / (2ai) , lim L^A^A = / 1 ( 2ai ), (104) 
p-5.00 p p^oo 2p oai 
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Table 1: Critical, T c , and self-dual, T* . temperatures. For p > 5, there are two critical 
temperatures. The lowest one, Tg^-p goes to zero when p — > oo, as Tg^ T ~ 1/p 2 , and the 



highest critical temperature 


T$ T ~ 0(1). 




V 


T c 






[J/k B ] 


[J/k B ] 


2 


2/1b(1 + v^) 


2/\n(l + V2) 


3 


3/(21n(l + y/3)) 


3/(21n(l + y/3)) 


4 


l/ln(l + V^) 


l/ln(l + V / 2) 


6 




l/(21n(2cos(f))) 


large p 




2vr/p 



and from the asymptotic expansion of the modified Bessel functions 7o,i> Eq. 
(27), one gets the relation between the transfer matrix and direct couplings 

1 4vr 2 

K 2 . (105) 



2a x {K 2 ) 

One can then use Eq. (103) to obtain a relation between the two critical 
temperatures T^j and T^ T ( Ji = J 2 = J) 

*bT« t = ^rir- (106) 

r ruB-L BKT 

The Peierls argument developed in Appendix E, on the other hand, provides 
a rigorous scaling for the lowest transition temperature, TgKT = 0(X/p 2 )i as 
p — > oo. Thus this lowest critical temperature vanishes for the classical XY 
model (where a broken-symmetry phase is not allowed) and, according to Eq. 
(106), TgKT ~ 0{\) for large p. The self-dual temperature T* has its own 
"intermediate" scaling with p, k B T* = (2ir/p) J, so that it also vanishes in 
the p — > oo limit. This is to be expected since the XY model is not self-dual, 
nor has a natural self-dual approximation. 

To understand what makes p > 5 different from p < 5 and explain the ap- 
pearance of the intermediate critical phase, one needs to analyze the nature 
of the topological excitations. For p > 5 there are two main types of topolog- 
ical excitations (see Fig. 5): (i) domain wall excitations that dominate the 
low-temperature physics, and (ii) discrete vortex-like excitations of relevance 
in the critical and high-temperature phases. The key distinction between the 
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Figure 5: Two types of topological excitations: domain wall (left panel) and discrete 
vortex-like (right panel) excitations. Integer numbers s = 0,1,2,3 indicate the angle 
variables 6 = 2ns/p. For p > 5 the energy cost of domain walls depends on | As| as opposed 
to the 2 < p < 4 case where the cost is independent of |As|. Vortex-like configurations 
start appearing at p = 5, and the example above shows a vortex in Z5 of strength k = 1. 



two is that domain walls exist for any p, while vortex-like excitations do not 
exist for 2 < p < 4, becoming manifest only for p > 5. Also, if 2 < p < 4 the 
energy cost to create a domain wall is independent of |As| = \s- m — s out |, with 
•Sin(out) indicating the angular configurations at the two sides of the wall. This 
changes for p > 5, allowing for twists of the spin of size |As| = 2, ■ ■ • , (p — 2). 

In the Peierls argument provided in Appendix E, the upper bound for 
the probability of having a domain wall corresponds to a change of orienta- 
tion between two domains of (±27r/p) (or equivalently |As| ~ 0{\)). Such 
a change always appears in all domain walls in systems with p = 2,3,4. 
This, however, is not the case for p > 5 where there exist general domain 
wall topologies that do not allow for a uniform twist of the angle between 
neighboring domains. In such instances, |As| can be 0(p) and, as shown in 
Fig. 5, two types of excitations are generally possible. More precisely, in 
p > 5 systems, a vorticity arises. The topological invariant characterizing 
these configurations, that we call discrete winding number k, is given by the 
circulation sum 




(107) 
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taken around an oriented loop T, with the argument A8 rr > G [— it, it) given 
by 

A8 rr > = (9 r , - 6 r ) mod(27r). (108) 

In contrast, the sum in Eq. (107) should not be taken mod(27r) but just as 
an ordinary sum of real numbers (otherwise it would vanish identically). To 
make this definition lucid, consider a (p — 5) configuration such as the one 
shown in Fig. 5, with = 2irs/p, and a loop T. Herein, the circulation sum 
explicitly reads (A0 O i+A0i 2 +A023+A0 3 o) = 27r/5+27r/5+27r/5+47r/5 = 2vr. 
Thus, the configuration in Fig. 5 has a vortex of strength k — 1 at its origin. 
As 0o — 03 = —6ir/5, we set A0 3O = 47r/5. This shift in the value of A0 rr r 
(that must lie in the interval [— 7r, 7r)) leads to the non-zero value of k in this 
case. 

We may now use energy-versus-entropy balance considerations to argue 
for the relevance of these topological excitations in establishing the two phase 
transitions. The Peierls argument presented in Appendix E rigorously es- 
tablishes that domain walls oriented relative to one another by the minimal 
energy cost (i.e., twists of (±2n/p)) are responsible for the existence of a 
low-temperature ferromagnetic broken Z p symmetry phase. Both the energy 
penalties and entropic costs associated with such minimal cost domain walls 
scale with I (the domain wall length), and the analysis leads to a transition 
temperature that behaves as ksT^j ~ (1 — cos ^y)J- The second transition 

temperature, T^ T (p > 5), is associated with the proliferation of vortex-like 
excitations, and as indicated in Appendix E should scale as Tg KT ~ 0{\). 
Note that this energy-versus-entropy balance argument does not rely on the 
existence or non-existence of the self-dual property of the model. 

It is important to mention that while the physics of the low-temperature 
phase is associated to the exact discrete Z p symmetry, the existence of vortex- 
like excitations is directly related to the emergence of the continuous sym- 
metry Ucf, unveiled in Section 3.5. This U(l) symmetry becomes more exact 
at high temperatures (T > T* or 2a\ > Ki), for a fixed p > 5, or it is exact 
at any temperature when p — > oo. Thus, the physical origin of the critical 
phase and the extended universality concept introduced in Ref. [19] is simply 
our emergent continuous symmetry. 

What is the nature of the phase transitions when p > 5? Given the 
current debates [22], it is important to say what we mean by a "BKT-type 
phase transition" . We simply mean a transition characterized by an essential 
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singularity in the free energy (or the ground state energy in a quantum 
model). This includes those cases where there is an essential singularity but, 
for instance, the correlation function exponent rj ^ 1/4 (1/4 is the exponent 
for the XY model [11]). It is very difficult to prove analytically the existence 
of an essential singularity but numerical simulations seem to indicate that for 
p > 5 the two phase transitions are continuous with continuous derivatives 
[22] , supporting the BKT scenario. Moreover, our new self-duality argument 
proves that the two transitions must be of the same nature [13]. In other 
words, if there is an essential singularity (the function and all its derivatives 
remain continuous) in the free energy at Tg^ T , then, there should be the 
same type of singularity at Tg^ T [13]. This, of course, does not mean that 
the self-duality fixes the value of, for instance, rj, to be the same at the two 
transition points. 

In what follows, we will show that when p > 5, the temperature region 
^bkt — ^ — ^bkt mus t be critical (massless) with algebraic correlations. The 
assumptions are: (1) The phase transitions at Tg^ T and T B ^ T are continuous, 
and (2) for large separations \r— r'\, the two-point correlation function G(\r — 
r'\,T) = (cos(6 r — 9 r ')) is of the canonical Ornstein-Zernike-like form 

e -\r-r'\/Z 

G(\r - r'\,T) ~ A ^ + M 2 , (109) 

7* — 7*' r ' 

with M representing the order parameter, i.e., magnetization, £ the correla- 
tion length, rj an anomalous exponent, and A an amplitude. In principle, all 
of these quantities (M,£,r), and A) are functions of temperature T. 

We start by demonstrating that G(\r — r'\,T) is a monotonically de- 
creasing function of T for any p. To this end, we derive in Appendix F a 
Griffiths'-type inequality as in general ferromagnetic systems [34]. Assume, 
for simplicity, the uniform case — K > 0. Then, it is straightforward to 
show (see Appendix F) that 

d K G(\r-r'\,T)>0, (110) 

or, equivalently, dxG(\r — r'\,T) < 0, which proves that G is a monotonically 
decreasing function of temperature (0 < G(\r — r'\,T) < 1). 

Now, if G(\r — r'\, T) is monotonically decreasing with temperature T for 
any fixed separation \r — r'\ then, in the absence of magnetization (i.e., when 
M = 0) the correlation length £ must be a monotonically decreasing function 
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of temperature. The proof of this assertion is trivial. Consider two rather 
general temperatures in this region T a , T b , such that T^ T < T a <T b < Tg^ T , 
then the ratio of the corresponding asymptotic correlation functions is 



because of the monotonicity property of G. This is only possible if < £ a 
for otherwise for large |r — r'\, the ratio of the two correlation functions 
would diverge exponentially in \r — r'\. As T a and T b were rather general 
temperatures, it follows that £(T) is a monotonically decreasing function of 
the temperature T so long as the magnetization M = (as it is for T > Tg^ T ). 
Therefore, a divergence of the correlation length at Tg^ T and Tg^ T implies 
that within the entire interval Tg^-p <T < T^ T , the correlator G is algebraic 
in |? — r'\. We thus proved that there must exist a power law, critical phase, 
between 2g^ T and Tg^ T . 
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Appendix A. Exponential of shift operators 

In this appendix we collect some useful formulas to compute exponentials 
and logarithms of shift operators. Let us start with the operator L z , the 
infinitesimal generator of translations on the circle. We have the general 
relation 



Jo 

Our goal is to compute a as a function of b, and vice versa. The first step is 
to notice that the Fourier transform operator 



G(\r-r'\,T b ) _ A b ^-r^-^) 
G(\r - r'\,T a ) ~ A a \r - r'^-va 



< 1, 



(111) 




(A.l) 




(A.2) 



puts L z (and the expression of Eq. (A.l)) in diagonal form, 




(A.3) 
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This, together with the orthogonality relation 2ti5(9' — 6) = Yl n & 
leads to 



b(0) 
a{B) 



— ' 



2w « 



2tt 



e ltfn In 



■2tt 



d9'b{e')e- Wn 



,\(9'-9)n 

) 

(A.4) 
(A.5) 



In Section 3.1 we introduced a diagonal matrix U, and a shift operator 
that describes translations in a p-points discretization of the circle. This 

shift operator plays a role similar to that of L z (actually, e~' 6Lz ). Now we 

have the general relation 



V p_1 a V r ' 



P-1 



h ^ 



(A.6) 



m=0 



that should be compared to Eq. (A.l). Our goal is to find closed-form 
expressions for the coefficients a m in terms of b m and vice versa. The unitary 
transformation that diagonalizes is now given by the discrete Fourier 
transform of Eq. (62). Putting these pieces together, we get the solution to 
our problem, 



1 

j m — — tr 
V 



Jjrn e ELo a ' Utl 



a 



1 

— — tr 

V 



p-i 

1=0 



• (A.7) 



As seen by expanding the trace, these equations are closely related to (A.4) 
and (A.5). 

In physical applications, the a m are Hermitian-symmetric, a p - m = a* m (to 
guarantee that Y7m=o a mV m is a Hermitian operator), and the b m are real 
and positive. Thus, it is convenient to assume that both set of coefficients 
satisfy a p _ m = a m , b p _ m = b m , and the relations between them simplify to 



1 ^ 

1 p ~ l 



COS 



COS 



s=0 



,2-Kms, 
V 

2ixms, 
V 



'p-1 



ME* 



;COS( 



.1=0 



,2ttIs. 

p 



(A.8) 
(A.9) 
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These are the expressions that are most useful in physical applications. 

Suppose next that b m = e Ku ( m \ where K is a positive constant, and u(m) 
is a real function of m — 0, • ■ ■ ,p — 1 (for example, u(m) = cos(^p) for the 
classical p-clock model). We would like to study the behaviour of the a m to 
next-to-leading order in K, in the limit that K grows very large (this could 
happen at low temperature). Notice that in this limit 

y bl C0S (— ) « e *«(°) ( 1 + 2e ^( 1 )-<°)) cos(— )) (A.10) 

to next-to-leading order, assuming that the inequalities 

> (u(l) - u(0)) > (u(2) - u(0)) > • • • (A.ll) 

hold. The factor two in Eq. (A.10) is due to the symmetry u(p — I) = u(l). 
Replacing expansion (A.10) into Eq. (A. 8) leads to 

a m w Ku(0)6 m , + e K ^~ u ^\6 mA + K -+ oo, (A.12) 

where we have used ln(l + x) ~ x, and the orthogonality relation 
p— l 

lv^ ,2'Kms. ,2irsl. 1 . . . . ,. 10 , 

- > cos( ) cos( ) = - [d m> i + Om,p-l) ■ (A. 13) 

p L — ' p pi 

s=0 



Appendix B. Duality of the XY model to q-deformed bosons 

In this section we study a duality that illustrates the essential differences 
between compact, 9, and non-compact, x, degrees of freedom. The algebraic 
tool of choice is the g-oscillator algebra [35], specified by a positive real 
number q, a creation operator a', its Hermitian conjugate a, and a Hermitian 
operator h, satisfying 

[n, a] = —a, [h, a'] = a , aa^ — ga^a = q~ n . (B-l) 

If g = 1, this algebra reduces to the standard harmonic oscillator algebra, that 
is isomorphic to the Heisenberg algebra of translations on the line, [x,p] = i, 
provided a — (x + \p)/ \/2, a) = (x — \p) / a/2. It was pointed out in Ref. [36] 
that the mapping 

L z ^ -n + ln v/2sinh(l), (B.2) 
e' § ^ y/2siuh{2)ae- fL , e~' § ^ v/2sinh(2) e"V, 
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affords a representation of the algebra of translations in the circle [L z , e ±id ] = 
±e , provided that we set q = e~ 2 in Eq. (B.l). 

Clearly, we can extend the mapping of Eq. (B.2) to a duality isomorphism 
$d f° r the XY model. The dual transfer operators read 

T? = y2sinh(l) I d6 e K2COS0 e ,9h \ (B.3) 
i=i 

N-l 

if = n ex p 

i=l 

To compute Z^y, one should take the trace in the eigenbasis of ?v This basis 
is described in Ref. [35]. 

This description of the XY model in terms of g-deformed bosons with q = 
e~ 2 suggests that the algebra of Eq. (B.l) affords a continuous interpolation 
between the XY model and ordinary phonons (characterized by q — 1), but 
this is not the case: The XY model belongs to a representation of the algebra 
of Eq. (B.l) that is inequivalent to that describing phonons (i.e., canonical 
bosons). The reason is that Eq. (7) is not enough to specify the algebra of 
translations in the circle. We must also have that 

The mapping of Eq. (B.2) will respect this constraint only if a, a) satisfy 

aa ] - q' l a)a = (B.5) 

at least for q = e -2 , including the relations listed in Eq. (B.l). But the 
resulting set of four relations becomes inconsistent at q — 1. This shows that 
the g-oscillator algebra cannot interpolate continuously between canonical 
bosons and compact excitations. 

Appendix C. The Villain and its dual solid-on-solid models 

The Villain model [11] 

{"(r, M )} i e r} 



K l smh{2){a l+l e-^ +1+hl) al + h.c) 



r u=l,2 



6L — 2irn 



(c.i; 
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was introduced in Ref. [16] to provide a Gaussian approximation to the XY 
model that preserves the essential property of compacticity, and is a good 
approximation at sufficiently low temperatures. We now show, by using our 
bond-algebraic approach, that it is dual to the solid-on-solid (SoS) model of 
the roughening transition, 



-Zsost-KjJ = E exp 

{m r } 



E E 

r fi=l,2 



(C.2) 



characterized by integer- valued degrees of freedom m r £ Z [11]. We work 
directly in the thermodynamic limit, N — > oo, to avoid dealing with boundary 
terms. 

The transfer operator T$ s — T2T1 for the SoS model can we written as 
Tl = ft e ^^ +1 -^) 2 , T 2 = n E ^ m2 i?! m , (C.3) 

j i m 

in terms of the operators X^ Ri, R\ defined in Eq. (32). Now, however, 
i £ Z labels the sites of an infinite straight line. The duality of bond algebras 

Xi-Xi-i A L z „ Ri -% eM+i-fc), 4 A e-Kft+i-^*) ( C .4) 

affords a dual representation of TsoS, 

if = H e^^S Tf = J] E e^V'^ 1 -*)™, (C.5) 

i i m 

in terms of compact degrees of freedom. The next step is to compute Z® oS = 
tr [(TlfT] )^] in the basis introduced in Eq. (9). 
is already diagonal in that basis 

T 2 D |#) = ] [ E e^ m \-' t{ - 0i +^- 6i ^ m \9). (C.6) 

i m 

At this point we could proceed by analogy to previous sections and rewrite 

"2 2 



m „—\0m 



this expression in terms of an interaction potential Vk(0) = — In e~ m e 
but this will not turn out be the most convenient approach. Instead, let us 
proceed to compute the matrix elements of I\ D . This task reduces to com- 
puting the matrix elements of a one-body operator, 

(O'fi^m = -Y (C.7) 

mi 
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which results from recalling that the orthonormal states of L z ^ are the plane 
waves (9i\rii) = e' 6irH j ' '. Notice that the function e% x2 is the Fourier 

transform of e^K j\/K. It then follows that we can use Poisson's summation 
formula to write 



E 



e T m i e~' 9rrii 




E 



■fit,- 



e 2A- 



(C.8) 



Putting all the pieces together, we obtain 

N 



'SoS 



%) En n e-p 

{e r } r [1=1,2 m 



^(9 r+e „ - 9 r - 2ixmf 



(C.9) 



2tt 
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S ex p 

{«(r,M)l IM 



EE^< 



9 r - 2nn {r ^ ) ) 2 



r ^=1,2 

The last expression is exactly (2ir / K2) N Zy [K^], and thus the Villain model is 
dual to the SoS model. Notice the reciprocal relation between the couplings: 
The Villain model is strongly coupled only if its dual SoS representation is 
weakly coupled. 



Appendix D. The p-clock model for p = 2, 3, and 4 

Let us start with the simplest p = 2 case. Then, Ui = U\ = of and 
Vi = V} = af, and the transfer matrix T p = T 2 Ti of Eq. (64) reduces to 



N-l 



N 



Tl =] e *-^+i, T 2 = J] {e K - + e- K *of) • 



(D.l) 



i=i 



i=l 



This finite Ising model is self-dual up to boundary corrections. The substitu- 
tion T\ — > e Kia NTi renders the model exactly self-dual for any N [13]. 
If p = 3, then V} 2 = Vi, and T 2 becomes 



N 

T 2 = n(^ 2 + e"^(^ + ^)). 
i=i 



(D.2) 



Ti is just as in Eq. (64), with Us appropriate for p = 3. It follows that if we 
introduce the boundary correction Ti — > e^^ UN+u ^Ti [13], then T 2 Ti 



and Ti 



T 2 , rendering Z p exactly self-dual for any N. 
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The case p = 4 is special because it can be mapped onto two decoupled 
Ising models [32]. Since V? = V} 2 in this case, T 2 reads 

N 

T 2 = Y[{e Ko +Vi + V? + e- K *V?). (D.3) 

i=l 

Moreover, it is easy to check that the operator 2Cj = (I+erf^) + (1 — of^crf^ 
(known as a controlled-NOT gate in quantum computation) maps 

(D.4) 

and thus it follows that C = Y\a=x C« maps 

N-l 

CTxC = JJ e^Ki^+i+^,i+i>, (D.5) 

i=\ 
N 

CT 2 C = JJ(e* + e -^^)(e* + e-^), (D.6) 

i=l 

that clearly defines two decoupled Ising models, with couplings that are half 
of those of the p = 2 model. In particular, the p = 4 clock model is exactly 
self-dual provided T x -> e ^Kiv+^,iv) t x . 

Appendix E. Peierls argument for the p-plock model 

We now use the Peierls argument to prove that there should be a broken 
symmetry phase (low-temperature ordered phase) in the p-clock model on 
the square lattice. The proof establishes the existence of a phase transition 
at a temperature T^> below which global Z p symmetry is broken. For large 
p» 1, T« = 0(l/p 2 ). 

Specifically, our objective is to show that if uniform boundary conditions 
pertaining to one of the clock states 9 = 2ns /p, with < s < p — la fixed 
integer, are applied on the boundary of the square lattice, then there prov- 
ably exists a temperature Tp e i er is > such that for temperatures T < Tp e i er | s , 
spontaneous symmetry breaking (SSB) of the global Z p symmetry arises. (In 
the context of our discussions thus far, Tp eier | s < T^; asymptotically, for 



tt i^o f a U ~ ,a li\r^ 

Ui = e'4Q^ JCi, 

Vt + V? = Q(a? i + <7*)Q, Vf = 
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large p, both temperatures scale as 1/p 2 .) By SSB in this context, we refer 
to the lifting of the symmetry triggered by applying the uniform boundary 
conditions at spatial infinity That is, when the aforementioned boundary 
conditions are introduced then, for T < Tp eier | s , the probability distribution 
V(8q) for the angular orientation of the spin at the origin So is not symmet- 
ric between the p possible values of 9q. In particular, we will demonstrate 
that V(9q) is maximal when 6q has an orientation that matches that on the 
boundary, 9^. In other words, for temperatures T < Tp e i er [ s , 

V(6 = O > -. (E.l) 
p 

To prove this inequality, we note that 

V(6 ^ 6^) < Prob(outer domain wall T), (E.2) 

where Prob(X) denotes the probability of the set of events X. The domain 
wall is defined as the boundary between differently oriented spins. The logic 
underlying Eq. (E.2) is clear: if 9 ^ 6^ then, by its very definition, at least 
one domain wall must separate the spin at the origin from the spins on the 
boundaries of the lattice. 

We now will bound the probability of having a particular domain wall. 
Specifically, let us denote by {C a } the set of configurations that have T as the 
outer-most domain wall surrounding the origin. That is, T separates spins 
with an orientation 8 = Boo from those having another (uniform) orientation 
6- m . [Note that, generally, more than one domain wall may be present and 
thus 9 m need not be the same as 6 .] The upper bound on the probabilities in 
Eq. (E.2) is a sum over the probabilities of having such different outer- most 
domain walls T. We furthermore define the partition function 

Z r = e M-f3E a ], (E.3) 

where E a = E(C a ) is the energy of the spin configuration C a . Thus, Z r is 
smaller than the total partition function Z p of the system. This is so as Z-p 
contains only a subset of the Boltzmann weights appearing in Z p . That is, 
in Eq. (E.3) we sum only over spin configurations with at least one domain 
wall surrounding the origin. 

We next define a new configuration C a formed by rotating all of the spins 
inside T by a uniform angle A8 such that the outermost domain wall that 
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surrounds the origin is removed. That is, for all spins S r that (i) lie inside 
the region bounded by the domain wall T, we perform the transformation 
9 r — > (9 r + A6) with an angle of rotation 

A# = 6>oo - 9- m = —As, (E.4) 
V 

where As is an integer, (ii) All spins lying outside the domain wall T have 
an orientation 9 = these spins are not rotated. When present, any other 
more internal domain walls will remain unchanged by this uniform rotation 
of all the spins inside V. In order to bound, from above, the probability of 
having an outermost domain wall T, we now consider {E a = E{C a )) 

Z f = Y^exp[-/3E a ]. (E.5) 

{<?«} 

The probability of having the domain wall T is fixed by the ratio of the 
sum of Boltzmann weights associated with having the domain wall T divided 
by the sum of Boltzmann weights associated with all spin configurations (i.e., 
the partition function Z p ). As Zf contains a sum only over a subset of all 
Boltzmann weights that appear in Z P1 we have 

Z v Z T e-^d + e- pEc 2 + ■■■ 

Probfouter domain wall V) = — — < — — = — . (E.6) 

V ' Z p ~ Z f e - fiE 0i + +... v ; 

The smallest energy difference between a configuration C a and C a is bounded 
by 

E Ca -Ec a >£(l-cos—), (E.7) 

where £ is the length of the domain wall T (exchange constants are set to 
unity, J M = 1, n = 1,2). As Eq. (E.7) applies to all configuration pairs C a 
and C a that appear in Eq. (E.6), 



e 

we have that 



Z T 



42 



It is important to emphasize that when the bound of Eq. (E.9) is saturated, 
| As | = 1. [Physically, for p ^> 1, only such domain walls (as opposed to far 
more energetically prohibitive domain walls with |As| = 0(p)) may appear 
at sufficiently low temperatures (T < 0(l/p 2 )).] 

Returning to the probability that (at least) one domain wall surrounds 
the origin in Eq. (E.2), we have that 

Prob(outer domain wall T) < N e D e , (E.10) 

t 

with Ni denoting an upper bound on the number of domain walls of perimeter 
i that enclose the origin and De an upper bound on the probability of having 
a domain wall of length £. Inserting Eq. (E.9) while taking note of an upper 
bound of 4 x 3 e ~ l on the number of non-backtracking walks of length i on 
the square lattice, and an upper bound of (^/4) 2 on the maximum number 
of initial starting points for a walk of length I that surrounds the origin, we 
have 

Prob(outer domain wall r) < ^ [(£/4) 2 x 4 x 3 ^V mi ~ cos f > 

= w(P,p) (E.ll) 

(the minimal domain wall on the square lattice has length £ = 4). The 
function w is trivially bounded by performing the summation over all natural 
numbers I 

wfap) < y ^! e -w-co S f) = x{i+x ) ay } (E 12) 

KH,FJ _ ^ 12(X-1) 3 V ' V ' 

■i-U j8(l-cos — ) /n 

with x = e p ' /S. 

In performing the summation in Eq. (E.12), we assumed a sufficiently 
low temperature so that x > 1, and a) is a monotonically decreasing function 
of (5. Notably, w can be made arbitrarily close to zero for large enough 
f3. Let us denote by /3p e \ er \ s the solution to the equation w;(/3p e i er | s ,p) = ^y-. 
Then, for > /3p e ieris, the probability of the spin at the origin being the 
same as that on the boundary is V(6o = 6^) > 1/p. In other words, for 
T < Tp e i er is, we clearly have SSB. It is important to emphasize that Tp e ieris 
is only a lower bound to the transition temperature, and the actual SSB 
occurs for > Tp. eier | s . An estimate for Tp eier | s resulting from this analysis 
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is Tp e j er |s ~ (1 — cos y )/ In 6 (~ 0(l/p 2 ) for large p). As in the lower bound 
derived herein, the energy cost for a domain wall is anticipated to determine 
the actual ordering temperature. This bound is rigorous. Physically, in Eq. 
(E.10), the logarithms of the two terms Nt and Di capture, respectively, 
bounds on the entropy and energy costs associated with domain walls of 
length £. 

Discrete vortices such as the one shown in Fig. 5 with a typical change of 
angle A8 = 0(1) (or |As| = 0(p)) across the intersecting domain walls that 
extend over a linear distance £ may entail, for all p > 5, an energy cost that 
scales as £. This is to be contrasted with the minimal energy penalty asso- 
ciated with a difference in angle of |A0| = In jp for which the corresponding 
energy penalty as £/p 2 (and that physically sets the bounds that we derived 
in the Peierls argument above). Thus, from energy- versus-entropy balance 
considerations, the temperature below which it is unfavorable to have vor- 
tices is Tvortex ~ C(l) ( or of order J): the energy for such domain walls scales 
as £ as does the entropy associated with a network of possible intersecting 
domain walls that have a total length t. 

Appendix F. Proof of monotonicity of the correlation function G 

We now prove Eq. (110) for large (yet finite) lattices. In finite size sys- 
tems (no matter how large), there are no thermodynamic phase transitions. 
Thus, the free energy and all its derivatives (including, in particular, the 
two-point correlation function G(\r — r'\,T)) are analytic for all values of /3. 
We will prove monotonicity by expanding (the analytic) G as a power se- 
ries in /3, which is everywhere convergent, and illustrate that the coefficients 
multiplying each power of /3 are non-negative. The uniformity of the sign of 
all contributions to the series coefficients follows from repeated applications 
of the identity 

p-i 

S p( n ) = ^2 e '~ S = P S n,0, (F.l) 
s=0 

with the Kronecker delta above defined mod(p). That is, 5 ,o = $±p,o — 
<5±2p,o = = otherwise it vanishes. 
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Longhand, the correlator G(\r — r'\,T) is given by (K = (3 J > 0) 



cos(0 r — 9 r r) exp 

(cos(0 r — 9 r r)) = ^—^ 



x fi=l,2 



Z P [K] 



(F.2) 



with 9 r = 2ns r /p. 

We Taylor expand the argument of the exponential, i.e, exp[/3A] = 1 + 
(5 A + (j3A) 2 /2\ + ■ • • , in both the numerator and the denominator of Eq. 
(F.2), and may represent pictorially the expansion terms by Feynman-type 
diagrams. In this scheme, each appearance of cos(9 x > — 9 X ) relates to an 
internal propagator linking sites x and x'. Similar to perturbative schemes 
in the continuum, where the number of Feynman-type "bubble diagrams" in 
Z V \K\ with nearest-neighbor links that share common vertices with any given 
"connected" diagram which contains both the sites r and r' is negligible, 
there is a near cancellation, in Eq. (F.2), of all "bubble diagrams". What 
remains from the ratio of Eq. (F.2) is the sum of all diagrams in which sites 
r and r' are linked to each other by the line stemming from cos(0 r — 9 r >) 
as well as internal lines which pass through the points x = 1, 2, • ■ ■ , m. Up 
to (inherently positive) symmetry factors, the numerical value of any such 
resulting diagram is given by a sum of the form 

p-i 

/ g*ri+tra+-+tia+-+t mr i cos (0 p _ p/ )( cos (0 r - 9±)) Ul 

S r ,Si,S2,— ,S m ,S r /=0 

x (cos(# r - 9 2 ))^ ■ ■ ■ (cob(0i - 9 2 )) tl2 ■ ■ ■ (cos(0 m - P O )*"**' , (F.3) 

where the integers tab > represent the number of lines linking sites a and 
b. It is simple to show that sums of the form of Eq. (F.3) are manifestly 
non-negative. Replacing cos(0 a — 9b) by (exp(i(0 a — 9 h )) + exp(i(0 fe — 9 a )))/2, 
a sum of exponentials with positive weights, Eq. (F.3) reduces to a sum of 
individual products with positive weights, each being of the form 

n > o> ( f - 4 ) 

x=r,l,2,--- ,m,r' 

with n a set by sums of the powers t a b in Eq. (F.3). Then, it follows that 
when the ratio of Eq. (F.2) is expanded in f3, i.e., G(\r — r'\,T) = Ylt a tP t i 
the prefactor a t multiplying each individual power /3* is non-negative. As 
such, G(\r — r'\,T) is manifestly monotonic in (3, i.e., d@G(\r — r'\,T) > 0. 



45 



References 

[1] R. Potts, Proc. Camb. Philos. Soc. 48, 106 (1952). 

[2] F. Y. Wu, Rev. Mod. Phys. 54, 235 (1982). 

[3] V. L. Berezinskii, Sov. Phys. JETP 32, 493 (1971). 

[4] J. Kosterlitz and D. Thouless, J. Phys. C 6, 1181 (1973). 

[5] J. Kosterlitz, J. Phys. C 7, 1046 (1974). 

[6] J. V. Jose, L. P. Kadanoff, S. Kirkpatrick, and D. R. Nelson, Phys. Rev. 
B 16, 1217 (1977). 

[7] P. Minnhagen, Rev. Mod. Phys. 59, 1001 (1987). 

[8] D. B. Abraham, Surface Structures and Phase Transitions - Exact Re- 
sults, in Phase Transitions and Critical Phenomena Vol. 10, edited by 
C. Domb and J. L. Lebowitz (Academic Press, 1986) p. 1. 

[9] D. R. Nelson and Halperin, Phys. Rev. B 19, 2457 (1979); A. P. Young, 
Phys. Rev. B 19, 1855 (1979). 

[10] S. Elitzur, R. B. Pearson, and J. Shigemitsu, Phys. Rev. D 19, 3698 
(1979). 

[11] H. Nishimori and G. Ortiz, Elements of Phase Transitions and Critical 
Phenomena (Oxford University Press, 2011). 

[12] E. Cobanera, G. Ortiz, and Z. Nussinov, Phys. Rev. Lett. 104, 020402 
(2010). 

[13] E. Cobanera, G. Ortiz, and Z. Nussinov, arXiv:1103.2776vl [cond- 
mat.stat-mech]. 

[14] Z. Nussinov and G. Ortiz, Phys. Rev. B 79, 214440 (2009). 
[15] Z. Nussinov and G. Ortiz, Europhysics Lett. 84, 36005 (2008). 
[16] J. Villain, J. Phys. (Paris) 36, 581 (1975). 
[17] R. Savit, Rev. Mod. Phys. 52, 453 (1980). 



46 



M. B. Einhorn, R. Savit, and E. Rabinovici, Nuc. Phys. B 170, 16 
(1980). 

C. M. Lapilli, R Pfeifer, and C. Wexler, Phys. Rev. Lett. 96, 140603 
(2006). 

J. L. Cardy, J. Phys. A: Math . Gen. 13 , 1507 (1980). 

J. Frohlich and T. Spencer, Comm. Math. Phys. 81, 527 (1981). 

O. Borisenko, G. Cortese, R. Fiore, M. Gravina, and A. Papa, Phys. 
Rev. E 83, 041120 (2011). This work presents a summary of current 
debates in its introduction. 

Y. Tomita and Y. Okabe, Phys. Rev. B 65, 184405 (2002). 

V. A. Malyshev and E. N. Petrova, J. Math. Sciences 21, 877 (1983). 

J. Schwinger, Quantum Mechanics: Symbolism of Atomic Measurements 
(Springer Verlag, Berlin, 2001). 

P. J. Davis, Circulant Matrices (John Wiley, New York, 1979). 

M. Henkel, Conformal invariance and critical phenomena (Springer Ver- 
lag, Berlin, 1999). 

Handbook of Mathematical Functions, edited by M. Abramowitz and I. 
A. Stegun (Dover, New York, 1972). 

V. B. Matveev, Inverse Probl. 17, 633 (2001). 

J. J. Rotman, An Introduction to the Theory of Groups (Springer Verlag, 
New York, 1999). See specially Chapter 11, Example 11.5 on page 347. 

C. D. Batista and G. Ortiz, Adv. in Phys. 53, 2 (2004). 

For p = 4 we may represent the four equidistant spin orientations on 
the unit circle in terms of the normalized two component spin vectors 
(xliXD/Vz at sites r with xl' 2 — ±1- With this representation, the 
isotropic partition function reads 



Z P=4 [K] = ^2 ex P 

{Xt 'Xt" } 



, (F.5) 



47 



which shows that the p = 4 clock model is equivalent to two decoupled 
Ising systems. 

[33] For p > 3, ground states that are also eigenstates of C , C |\&±) = ±|*±), 
can be constructed: 

(i) p G odd: 

K) = » 

* s ± ) = ± M' S ))/V2 , for 1 < s < (p - l)/2, 

(ii) p G even: 

*°+> = « , l*f ) = l< 2 > 

W£) = (|*5) ± l^r 5 ))/^ , for 1 < s < p/2 - 1 

[34] R. B. Griffiths, J. Math. Phys. 8, 478 (1967). 
[35] G. Rideau, Lett. Math. Phys. 24, 147 (1992). 

[36] K. Kowalski, J. Rembielinski, and L. C. Papaloucas, J. Phys. A: Math. 
Gen. 29, 4149 (1996). 



48 



